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We present predictions for the pseudo-rapidity density of charged particles produced in central 
Pb-Pb collisions at the LHC. Particle production in such collisions is calculated in the framework of 
fct-factorization. The nuclear unintegrated gluon distributions at LHC energies are determined from 
numerical solutions of the Balitsky-Kovchegov equation including recently calculated running cou- 
pling corrections. The initial conditions for the evolution are fixed by fitting RHIC data at collision 
energies ^/Jnn = 130 and 200 GeV per nucleon. We obtain dN^^~^^ / drj (y^sjvjv ~ 5.5 TeV)|,,=o ~ 
1290 ^ 1480. 



It has been suggested that the nucleus-nucleus colli- 
sions performed at the Relativistic Heavy Ion Collider 
(RHIC) at the highest collision energies of 130 and 200 
GeV per nucleon probe the Color Glass Condensate [ij 
regime of QCD governed by non-linear coherent phenom- 
ena and gluon saturation. This claim is supported by the 
success of saturation models 0, d, 0, [E [B| in the descrip- 
tion of the energy, rapidity and centrality dependence 
of the particle multiplicities experimentally measured in 
d-Au and Au-Au collisions. With collision energies of 
up to 5.5 TeV, the upcoming program in lead-lead colli- 
sions at the CERN Large Hadron Collider (LHC) is ex- 
pected to provide confirmation for the tentative conclu- 
sions reached at RHIC and to discriminate between the 
different physical mechanisms proposed to explain par- 
ticle production in high energy nuclear reactions (for a 
review of alternative approaches see, e.g., 

The phenomenological models in [3, 0, 0, 0] rely on 
the assumption that the saturation scale QsA that gov- 
erns the onset of non-linear effects in the wave func- 
tion of the colliding nuclei is perturbatively large 1 
GeV at the highest RHIC energies. Next, gluon pro- 
duction is calculated via the convolution of the nuclear 
unintegrated gluon distributions (ugd's) according to fct- 
factorization [8| . Under the additional assumption of lo- 
cal parton-hadron duality, the multiplicity in A-A colli- 
sions at central rapidity rises proportional to the satura- 
tion scale, dN / dri\ri=o (x Q^^, [9|. On the other hand, 
the growth of the saturation scale with increasing en- 
ergy (equivalently, decreasing Bjorken-a:) is determined 
by the perturbative BK-JIMWLK non-linear evolution 
equations [l^, [llj (for a complete set of references see 
[ij), thereby establishing a direct link between the ini- 
tial state gluon saturation dynamics and the experimen- 
tally measured hadron yields. The energy dependence of 
the saturation scale yielded by the BK-JIMWLK equa- 
tions at the degree of accuracy of their original derivation, 
leading-logarithmic (LL) in Q!sln(l/x) with a,, fixed, is 
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This 



A - 0.2 0.3 0, i, S, [13. Such deficiency of the theory 
has been circumvented so far by leaving A as a free pa- 
rameter, often adjusted to the empirical value A « 0.288 
obtained in fits to small-x HERA data in deeply inelastic 
lepton-proton scattering in the framework of saturation 
models 13l, Il5|. 

Higher order corrections to the BK-JIMWLK equa- 
tions have been calculated recently via all orders re- 
summation oiagNf contributions [l6|. Such corrections 
bring substantial modifications to the LL kernel, includ- 
ing running coupling effects, and result in a significant 
slowdown in the speed of the evolution, among other 
quantitatively important dynamical effects 17|. 

In this work we demonstrate that this partial improve- 
ment is sufficient to describe the energy and rapidity de- 
pendence of the multiplicities in Au-Au collisions at the 
highest RHIC energies. Then we extrapolate to LHC 
energies and present predictions for Pb-Pb collisions at 
y'SNN = 5.5 TeV. 

We start by solving the non-linear small-x evolu- 
tion equation for the dipole-nucleus scattering matrix, 
S'(y,r), including running coupling corrections 1^ 17 1: 



dS{Y, r) 
dY 



n[s]~s [s] 



(1) 



growth is too fast to be reconciled with the energy depen- 
dence observed in RHIC multiplicity data, which indicate 



where r is the dipole size and Y — ln(a;o/a;). The first, 
running coupling, term of the evolution kernel, 7?.[iS'], re- 
casts the higher order corrections that amount to a mod- 
ification of the LL small- a; gluon emission kernel, leaving 
the interaction structure of the LL equation untouched, 
whereas the second, subtraction, term, S[S], accounts for 
the new interaction channels opened up by the higher 
order corrections. Explicit expressions for both terms as 
well as a detailed explanation of the numerical method 
used to solve Eq. ^ are given in 17]] . The initial condi- 
tions for the evolution are taken from the semi-classical 
McLerran-Venugopalan (MV) model [HI, aimed at de- 
scribing the gluon distributions of large nuclei at mod- 
erate values of Bjorken-x, prior to the onset of quantum 
corrections. Thus, the initial dipole scattering amplitude. 
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J\f=l~S, reads 



JV{Y = Q,r) = l-cxp<^- 



(2) 



where Qq is the initial saturation scale. The constant 
e under the logarithm acts as an infrared regulator and 
A = 0.2 GeV. 
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FIG. 1: A: 
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for Qo = 0.5, 0.75, 1 and 1.25 GeV 



(from top to bottom), and for k = 0.5 (left) and k = e ^ 
(right). 

The speed of evolution, X = d\nQ^{Y)/dY, extracted 
from numerical solutions of Eq. ^ corresponding to dif- 
ferent initial conditions (Qo = 0.5, 0.75, 1 and 1.25 GeV) 
is plotted in Fig. [TJ For Y>0 the saturation scale is de- 
termined by the condition J\f{Y, r = l/QgiY)) — k, with 
K = 0.5 (left plot) and K = e~^ (right plot). These results 
show two remarkable features of the solutions. 

First, the running coupling corrections render the en- 
ergy dependence of the saturation scale compatible with 
the one indicated by the analysis of experimental data. 
Thus, the A values in Fig. [1] are slightly smaller than the 
one extracted from fits to HERA data, A = 0.288 (except, 
perhaps, for Qo ^ 0.5 GeV at small rapidities). On av- 
erage, they are compatible with A = 0.2 reported in as 
the optimal value to reproduce the energy and rapidity 
dependence of the multiplicities in Au-Au collisions at 
the highest RHIC energies. Second, they reveal the exis- 
tence of two very distinct kinematical regimes: At small 
pre- asymptotic rapidities the evolution is strongly depen- 
dent on the initial conditions. In particular, denser sys- 
tems, i.e. those associated to larger values of Qo, evolve 
more slowly due to the relative enhancement of non-linear 
effects with respect to more dilute systems. Such depen- 
dence on the nature of the evolved system is completely 
washed out by the evolution and, at high enough rapidi- 
ties, all the solutions reach a common speed of evolution. 
The onset of this universal scaling regime is reflected in 
Fig.[T]by the convergence of all the individual trajectories 
into a single curve for F > 15. The studies of more exclu- 
sive properties of the solutions carried out in 13, 17 1 sug- 
gest that the full scaling regime is reached at even larger 



rapidities, Y > 80. Moreover, sizable scaling violations 
have been detected in HERA data [1911 and in particle 
spectra in d-Au collisions at RHIG [20|. These observa- 
tions rise the question of whether the scaling ansatz that 
connects HERA and RHIC phenomenology through the 
universality property of the solutions is an adequate one 
at presently available energies. 

In analogy to [1, 0] , we calculate the pseudo-rapidity 
density of charged particles produced in nucleus- nucleus 
collisions within the fcf-factorization framework via: 
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where pt and y are the transverse momentum and ra- 
pidity of the produced particle, xi^2 = {Pt/ \/s) , Q = 
0.5 max{|pt ± fc* |} and b the impact parameter of the col- 
lision. The lack of impact parameter integration in this 
calculation and the gluon to charged hadron ratio are 
accounted for by the constant C, which sets the nor- 
malization. The nuclear unintegrated gluon distribution 
entering Eq. ([3]) is related to the inclusive gluon distri- 
bution, tp{x,k) (X ' ^^"^ given in terms of the 
dipole scattering amplitude evolved according to Eq. ([1]): 
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The relation between the evolution variable in Eq. ([T]) 
and Feynman-x of the produced particle is taken to be 
Y = ln(0.05/a;i^2) + AYg^. Since the relevant values of 
Bjorken-x probed at mid-rapidities and ^JSNN = 130 
GeV at RHIG are estimated to be - 0.1 0.01, the 
free parameter ^Yev controls the extent of evolution 
undergone by the nuclear gluon densities resulting of 
Eq. Il]) prior to comparison with RHIC data. Similar 
to 4], large- a; effects have been modelled by replacing 
tp{x,k) ip{x,k){l — x)'^. The running of the strong 
coupling, evaluated according to the one loop QGD ex- 
pression, is regularized in the infrared by freezing it to 
a constant value ajr — ^ at small momenta. Finally, in 
order to compare Eq. ^ with experimental data it is nec- 
essary to correct the difference between rapidity, y, and 
the experimentally measured pseudo-rapidity, ry. This is 
achieved by introducing an average hadron mass, m. The 
variable transformation, y(rj,pt,m), and its correspond- 
ing Jacobian are given by Eqs. (25-26) in 0]. Corrections 
to the kinematics due to the hadron mass are also consid- 



ered by replacing pf 



in the evaluation of 



xi 2. Remarkably, the optimal value found in comparison 
with data, m ~ 0.25 GeV, see Fig.[2]A, is in good quanti- 
tative agreement with the hadrochemical composition of 
particle production at RHIC. 

With this set up we find a remarkably good agreement 
with the pseudo-rapidity densities of charged particles 
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FIG. 2: Pseudo-rapidity density of charged particles pro- 
duced in Au-Au 0-6% central collisions at ^/snn = 130 GeV. 
Data taken from [ij. The solid lines correspond to Qo = 1 
GeV, m = 0.25 GeV, AY = 1 and 2:1,2 = (mt/^s) e*". The 
modifications to this central value considered are: Upper-left: 
m — 0.5 GeV (dashed line) and m — GeV (dashed-dotted 
line). Upper-right: AY = 3 (dashed line) and AY — 10 
(dashed-dotted line). Lower-left: Qo = 0.7 GeV, (dashed 
line) and Qo = 1.25 GeV, (dashed-dotted line). Lower-right: 
a;i,2 = (pt /y/s) e*^ (dashed line) . 



measured in 0—6% central Au-|-Au collisions at collision 
energies ^ysNN — 130 and 200 GeV. The comparison with 
data [l3| , shown in Fig. [2l constrains the free parameters 
of the calculation to the ranges: Qo ~ 0.75 1.25 GeV, 
m ~ 0.25 GeV and 3 > AYe^ > 0.5. These ranges deter- 
mine the uncertainty bands of the LHC extrapolation in 
Fig.H The best fits (solid hues in Figs. (HD and ([3])) are 
obtained with Qo = l GeV, m = 0.25 GeV and AY^v^l- 
The normalization constant, C, fixed at y/sJfN = 130 
GeV and f] — Q, is of order one in all cases. The line 
of argument that leads to these values is the following: 
First, the energy extrapolation from 130 to 200 GeV at 
central rapidities demands a moderate evolution speed 
A~0.2 0. From Fig. [1] that condition is met by either 
initial saturation scales Qo ~ 1 GeV and small evolution 
rapidities AY^^ ^ 3 or at asymptotically large rapidi- 
ties, AYf^y ~ 50, which are kinematically excluded. In 
the physically accessible range, the solutions close to the 
scaling region, i.e. for AY^v ~ 10, result in too narrow 
pseudo-rapidity distributions independently of the value 
of Qo, see Fig. [2j3. In the pre-asymptotic regime at fixed 
AYev ^ 3, those solutions corresponding to a Qo ^ 0.75 
GeV yield exceedingly broad distributions (see Fig.[2P). 
Thus, the energy and the pseudo-rapidity dependence in- 
dependently constrain the parameters of the gluon distri- 
butions probed at RHIC to the same ranges. This pro- 
vides the baseline for further evolution to LHC energies. 
In summary, these results indicate that the nuclear gluon 
densities probed at RHIC are in the pre-asymptotic stage 



of the evolution. This, together with the large values of 
the initial saturation scale required by data suggests that 
the saturation of gold nuclei at RHIC energies is not dy- 
namically generated by the evolution but, most likely, it 
is attributable to the nuclear enhancement factor that 
lies at the basis of the MV model, i.e., to the fact that 
the number of gluons in the nuclear wave function is large 
even at moderate energies due to the spatial superposi- 
tion of a large number of nucleon's gluon fields. 
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FIG. 3: Pseudo-rapidity density of charged particles produced 
in Au-Au 0-6% central collisions at -^snn = 130 and 200 
GeV and for Pb-Pb central collisions at ^sjvjv = 5.5 TeV. 
Data taken from 141 . The upper, central (solid lines) and 
lower limits of the theoretical uncertainty band correspond 
to (Qo = 1 GeV, AY = 1), (Qo = 0.75 GeV, AV = 3) and 
(Qo = 1.25 GeV, AY = 0.5) respectively, with m = 0.25 GeV 
in all cases. 

The extrapolation to LHC energies, done neglecting 
the differences between lead and gold nuclei and pre- 
sented in Fig. [21 is now straightforward and completely 
driven by the non-linear dynamics of gluon densities. For 
central Pb-Pb collisions we get 

Jl^^ ( 5.5 TeV, 77 = 0) - 1290 - 1480 , (5) 

with a central value corresponding to the best fits to 
RHIC data ~ 1390. These values are significantly 
smaller than those of other saturation based calculations 
H H, El, ~ 1700 2500, and compatible with the ones 
based on studies of the fragmentation region Such 
reduction is due to the lower speed of evolution yielded 
by Eq. ([T]) and to the proper treatment of pre-asymptotic 
effects, thereby going beyond the scaling ansatz. Impor- 
tantly, the prediction for the midrapidity multiplicity in 
Eq. (O is very robust against changes in the description 
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FIG. 4: Charged particle multiplicity in central Au-Au colli- 
sions at 77 = versus collision energy. Upper plot: Results ob- 
tained with the set up leading to Eq. ([5]) (band) and several 
modifications of it (see text). Lower plot: Power-law, a-^s'', 
and logarithmic, a+blns, fits to RHIC data at y^ijvjv ~ 19.2, 
64.2, 130 and 200 GeV. 

of particle production and the implementation of large- 
X effects. This is illustrated in Fig. |4lA, where the fol- 
lowing modifications to our set up have been considered 
(the itemization here follows the labeling in Fig. [JJ'V): 
a) Replacement of the ugd's in Eq. ^ by the modified 
gluon distributions h{Y, k) = fc^ </3(F, fc), as advocated 
in [2^. b) Regularization of the strong coupling at the 
value a/r. = 0.5, c) Removal of the (1 — a;)** corrections to 
the ugd's and d) Putting m — Q. The results obtained 
with these alternative configurations do not deviate from 
the uncertainty band given in Eq. ([5]), confirming that 
our predictions are mostly driven by the properties of 
small-x dynamics. Oppositely, our predictive power at 
large pseudo-rapidities, |ry| ^ 6, is lessened by the sen- 
sitivity of the evolution to the initial conditions and by 
our relatively crude implementation of large- a; effects (see 
Fig. [2P) , which are dominant in that region. 

Purely empirical parametrizations of multiplicity data 
of a large variety of colliding systems allow a logarithmic 
dependence on collision energy (see e.g. jlSl)- shown 
in Fig. HP, RHIC data by themselves do not differenti- 
ate between this and other functional forms like power- 
laws, negating any possibility to usefully constrain the 
expectations for LHC energies without further theoreti- 
cal guidance. Our results, similar to other calculations 
based on perturbative QCD, exhibit a power-law growth 
of the midrapidity multiplicity with increasing collision 



energy. The higher order corrections utilized here for the 
first time provide a richer physics input and result in a no- 
ticeably smaller power than previously estimated. This 
fact is crucial to obtain a good description of both the 
energy and pseudo-rapidity dependence of existing data 
and is the key ingredient in the extrapolation to higher 
energies. 
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